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Testing the Leutwyler-Smilga prediction regarding the global topological 
charge distribution on the lattice* 

S. Diirr ^ t 

'^Paul Schcrrcr Institut, Theory Group, 5232 Villigen PSI, Switzerland 

I give a sketch of my recent attempt to test the prediction by Leutwyler and Smilga according to which, for 
QCD in a finite box with Nf >2, the combination x = VEm indicates whether the net topological charge of the 
gauge background proves relevant (a; <§; 1) or irrelevant (a: ;^> 1) for physical observables. 



neously — provides a rich and fascinating subject: 
One expects to see both symmetry restoration 
phenomena and the onset of spontaneous sym- 
metry breaking (e.g. of the global SU{Nf)A sym- 
metry in QCD), if the box volume is taken suf- 
ficiently "small" or "large" , respectively. An ob- 
vious question is what sets the scale, i.e. by which 
standards does the box have to be small or large 
to trigger one or the other type of phenomena ? 
The naive guess is of course the lightest particle 
at hand, i.e. the mass M^^ of the Goldstone boson 
produced in the infinite volume limit. From this 
one would expect symmetry restoration pheno- 
mena to be pronounced for Af^L^l and SSB to 
become manifest for M-^L^^l. 

Another place where finite and infinite volume 
aspects of the theory meet is the question whether 
the total topological charge of the gauge field 
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Figure 1. Distribution of the naive topological 
charge I'nai (see footnote 1) in an ensemble of 
S'C/(3)-configurations generated with ^wiison at 
/3 = 6.1 (top) and in an ensemble of (full) QCD 
configurations (bottom) . Figures taken from [n| . 

1. INTRODUCTION 

QCD in a finite box — or, more generally, the 
finite-volume version of any field theory which 
breaks a continuous global symmetry sponta- 



*Talk given at QCD'OO, Montpellier, 
TE-mail: stephan.duerrOpsi.ch 



-12 July 2000. 



(which is a finite volume concept) proves relevant 
in physical observables (after extrapolating to in- 
finite volume) or not. From a lattice perspec- 
tive this question deserves attention, because of 
a technical subtlety with numerical QCD studies: 
In quenched simulations, excellent ergodicity of 
the sample w.r.t. the topological charge is easy 
to achieve (see top part in Fig. 1), i.e. the overall 
distributionFl is typically symmetric and centered 

^On the lattice a naive implementation of (|l|) — the re- 
sult being called u^^i — does not yield an integer, but a 
distribution which gets closer to a set of delta peaks near 
integer numbers only in the continuum limit a— >0. In the 
following, our interest is in the overall distribution, i.e. in 
the embedding curve of the distribution shown in Fig. 1 



about zero. In a dynamical run (one that includes 
the fermion determinant) the resulting distribu- 
tion may look like the one shown in the lower 
part of Fig. 1 — the simulation seems far from 
having achieved a good sampling of the topologi- 
cal charge. There are two reasons responsible for 
this: First, standard full QCD algorithms have a 
tendency to "get stuck" in a particular topologi- 
cal sector, if the quark mass is taken sufficiently 
light. This holds true both with the staggered |^ 
and with a Wilson type formulation H of the 
(dynamical) fermions. Second, because of the 
tremendous increase in costs (in terms of CPU 
time) per configuration, the total sample weight 
is typically smaller than in the quenched case. 
Hence the physical question is whether the full 
QCD sample shown in the lower part of Fig. 1 re- 
sults in unbiased measurements of Mjr, Ft^ etc. or 
whether the supposed overrepresentation of the 
sector with I' — S (see Fig. 1) is likely to affect the 
extraction of physical quantities. 

2. RESULTS OF THE ANALYSIS BY 
LEUTWYLER AND SMILGA 

A key point in the analysis by Leutwyler and 
Smilga [Q is that the topics raised in the intro- 
duction are closely interwoven and that the dis- 
cussion involves yet another parameterg 



ysm 



(2) 



which really decides whether the box is "small" 
or "large" . This new parameter is, in principle, 
independent of the naive parameter M^_L: In the 
regimes x 3> 1 and x ~ 1 with pronounced or mild 
SSB the box may be large [Mt^L ^ 1), interme- 
diate (Af^L ~ 1) or small {M^^L < 1) w.r.t. the 
conventional classification, whereas in the sym- 
metry restoration regime (x ^ 1) such a distinc- 
tion does not make sense anyways, because there 
the pion is not a useful degree of freedom. Note 
that in a lattice context the Leutwyler-Smilga 
(LS) classification refers to the mass of the sea- 
quarks (i.e. the dynamical quarks which influence 

^Here and below V denotes the four-volume of the box, 
S = — lim^^O lim\^— >oo('0V') (this order) is the chiral con- 
densate in the chiral limit and m is the (degenerate) 
sea-quark mass; both E and m are scheme- and scale- 
dependent, but the combination is an RG-invariant. 



the functional weight of a gauge configuration 
through the determinant), whereas the conven- 
tional one refers to the mass of the current-quarks 
(i.e. those from which the observable is built). 

The main result of the analysis by Leutwyler 
and Smilga Q is that there are two regimes of 
quark masses and box volumes for which the 
path-integral may be evaluated analytically and 
statements regarding the relevance of the topolo- 
gical charge v can be formulated. They find that 
in the symmetry restoration regime (x^ 1) where 
the fundamental description in terms of quarks 
and gluons is appropriate the partition function 
Z is completely dominated by the contribution Zq 
from the topologically trivial sector. In the oppo- 
site regime (x':S>l) where SSB becomes manifest 
(though formally the box volume is still finite) an 
effective description in terms of Goldstone degrees 
of freedom proves useful, but here two subcases 
must be treated separately: For M^L :i> 1 (i.e. if 
the pion stays inside the box) standard chiral per- 
turbation theory may be used with the obvious 
replacement ^jq^ -^ Y. p^+M^ ' i-*5- *^^ propa- 
gators should account for the "mirror copies". 
Since this perturbative evaluation assumes a sin- 
gle vacuum (sc. the one with u = 0) no statement 
regarding the z^-dependence can be made for this 
subcase. For M^L <C 1 (i.e. if the pion overlaps 
the box) the zero-modes provide the dominant 
contribution to the path-integral in the effective 
description ^, and Leutwyler and Smilga end 
up finding that the partition function is approxi- 
mately independent of v or, more precisely, that 
the distribution is [Q (a;>l,M^L<Cl) 

e ^<T^ with {i^^) = -^^^ . (3) 
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This is nice, but from a lattice perspective the 
statement about the z/-independence of the par- 
tition function in the large x regime is exactly in 
the wrong subcase — in numerical simulations the 
pion would always fit into the box. Hence, besides 
testing the LS-prediction for the small x regime, 
an obvious goal was to study the (/-dependence 
of both the partition function and selected ob- 
servables in the regimes with mild (x~l) or pro- 
nounced (x^l) SSB after reversing the "overlap- 
condition" into M^i>l 0. 
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Figure 2. Distribution of vah (top) and Vg^o (bottom) in the small (8x4 lattice), intermediate (18x6) 
and large (40x10) Leutwyler-Smilga regimes, respectively (QED(2) data for /3 = 3.4, r7T. = 0.09, 7V/ = 2). 



3. CHECKING THE NET TOPOLOGI- 
CAL CHARGE DISTRIBUTION 

In order to check the LS-prediction regarding 
the distribution of the net topological charge v in 
the regimes x-^l and x ^ 1 , one needs dynamical 
(i.e. unquenched) configurations. Given the alter- 
native to leave the subject to the phenomenologi- 
cal lattice groups which are able to generate full 
QCD ensembles or to study it in a suitably cho- 
sen model, I have opted for the latter possibility. 
For reasons to be suppressed due to limitations of 
space the massive multiflavour Schwinger model 
(i.e. QED in 2 dimensions with A^ > 2) is sup- 
posed to reproduce all qualitative aspects of the 
Leutwyler Smilga issue — on this point the reader 
is referred to |6|. My implementation uses the 
Wilson gauge action S'gauge = /3X](1~ ^os^n) and 
a pair of (dynamical) staggered fermions. 

The idea is to compare the three regimes x <C 1 , 
a; ~ 1, a: ^ 1 to each other using three dedicated 
simulations: Working at fixed /3 — l/{ag)'^ — 3.4 
and fixed staggered quark mass m= 0.09 (in lat- 
tice units), the three regimes are represented by 
the three volumes F = 8 x 4, 18 x 6, 40 x 10, respec- 
tively. From this setup the LS-parameter takes 
the values x ci 0.33,1.12,4.16, respectively, and 



the pion (pseudo-scalar iso-triplet) has a (com- 
mon) mass Mjr = 0.329 and hence a correlation 
length ^7r = 3.04 as to fit into the box. 

For the type of investigation I am aiming at 
configurations must be assigned an index v. I 
have implemented both the geometric definition 
i^gco = 27r X] I'^S ^n ^^'^ ^^^ field-theoretic version 
^fth = K i-'nai, ^'iiai = X] ^^^ ^^ with the rcnormali- 
zation k~1/(1 — (S'gaugo)//?^) 0- A configura- 
tion is assigned an index only if the geometric 
and the field-theoretic definition, after rounding 
to the nearest integer, agree. This turned out 
to be the case at a 99.9%, 98.8%, 88.2% level on 
the small/intermediate/large lattice, respectively. 
This fraction being so high means that in prac- 
tice an assignment can be done without cooling 
for the overwhelming majority of configurations; 
the remaining ones are just not assigned an index. 

As one can see from Fig. 2, the charge distri- 
butions found in the three runs seem to follow a 
general pattern consistent with the predictions by 
Leutwyler and Smilga: For x^l the distribution 
(an hence the partition function) is very much 
dominated by the contribution from the topologi- 
cally trivial sector, whereas for x S> 1 the distri- 
bution gets broad and seems compatible with the 
gaussian form (^) with {i''^)=x/2. 
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Figure 3. Scatter plots of Sformion^ — log{det{.^ + m)) (i.e. per continuum flavour) versus 5'gaugo (both 
normalized per plaquette □) on the complete sample (left) and for the sectors with i/ = and ^ — ±2 in 
the simulation covering the intermediate regime — each dot represents a configuration. 
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Figure 4. S'fcrmion versus iSgaugc together with the best linear fits for the reasonably populated topolo- 
gical sectors in the small (8x4 lattice), intermediate (18x6) and large (40x10) Leutwyler-Smilga regimes, 
respectively (QED(2) data for /3 = 3.4, m = 0.09, 7V/ = 2). 



4. ROLE OF THE FUNCTIONAL DE- 
TERMINANT 

Since the LS-issue is peculiar to the full (un- 
quenched) theory, an attempt to understand by 
which mechanism the three regimes differ from 
each other may lead one to investigate how the 
functional determinant or specifically its contri- 
bution to the total action per continuum-flavour 



5'fermion = - log(dct (.^ + m)) + COUSt 



(4) 



relates to the contribution from the gauge field 
and, in addition, how this might depend on the 
topological charge of the background. The idea 
is thus to study such a relationship sectorally, 
i.e. after the complete sample has been separated 
into subsamples with a fixed topological charge 
(or fixed \i^\), in the spirit of Ref. |8|. 

The first plot in Fig. 3 (which is for x ~ 1) 
tells us that there is a positive but weak correla- 
tion between S'fcrmion and 5*, 



gauge • 



From this one 
concludes that the functional determinant acts - 
roughly - like an effective renormalization of /3 



with a factor bigger than 1 . A key observation is 
that the correlation improves, if one separates the 
sample into subsamples with fixed \iy\, as is done 
in the r.h.s. of Fig. 3. Most notably the two sec- 
tors shown (i/ = and z^ = ±2) differ both in offset 
and slope of the best linear fit to 5'fermion versus 
5. 



gauge • 



This means that the functional determi- 
nant brings - in general - an overall suppression 
of higher topological sectors w.r.t. lower ones and 
a sectorally different renormalization of (3. 

Fig. 4 allows one to asses the strength of these 
effects as x varies: In the small LS-regime both 
the sectoral dependence of the renormalization 
factor and the offset between neighboring sectors 
are huge, i.e. for a;— >0 the functional determinant 
acts as a constraint to the topologically trivial 
sector. For x~\ both offset and sectoral depen- 
dence of the slope are very much reduced. In the 
large LS-regime there is only a minor overall sup- 
pression of higher topological sectors w.r.t lower 
ones, but the renormalization of (3 (and hence, 
in 4 dimensions, the lattice spacing in physical 
units) seems to be uniform for all sectors. 
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Figure 5. Scatter plots of the sectoral expectation values {L)±u of the Polyakov loop for the smallest 
and largest (reasonably populated) \iy\ in the small (8x4 lattice), intermediate (18x6) and large (40x10) 
Leutwyler-Smilga regimes, respectively (QED(2) data for /? = 3.4,to = 0.09, A(/ = 2). 



5. SECTORAL HEAVY QUARK FREE 
ENERGIES 

A question which has not been addressed so 
far is whether the statements by Leutwyler and 
Smilga regarding the i^-dependence of the par- 
tition function (which we found well obeyed) 
would carry on to observables. In other words: 
The question is whether in the small x regime a 
typical^ observable would be dominated by the 
contribution from the topologically trivial sector 
and likewise, whether the same observable would 
prove approximately independent of the topolo- 
gical charge of the background, if the analysis is 
performed in the large x regime. 



^Here we shall exclude observables that obey a topological 
selection rule, e.g. those that receive only contributions 
from the sectors with f = it 1 , as is the case for the fermion 
condensate in the massless limit. 



A quantity which is easy to determine in lat- 
tice studies and which might help elucidating the 
physical meaning of the LS-classification is the 
Polyakov loop, i.e. the trace of a chain of link vari- 
ables which winds once around the torus in the 
euclidean timelike direction. Its logarithm repre- 
sents (up to a factor) the free energy of an ex- 
ternal ( "heavy" ) quark which is brought into the 
system without possibility to influence it through 
back-reactions. As in the previous section the 
analysis shall be done sectorally, i.e. on classes 
of configurations with a fixed value of \i>\. This 
means that the physical (unseparated) expecta- 
tion value is understood as a weighted mean 



kl>o 



(5) 



of sectoral expectation values (i)|iy|, where the 



factor p\^\ reflects the combined weight of the sec- 
tors ±i/ in the appropriate histogram in Fig. 2. 
On a practical level, the separation (^ is facilita- 
ted by the fact that the ensemble is already repre- 
sentative in the sense of the full theory: In order 
to compute the ensemble average [i.e. the practi- 
cal version of the l.h.s. of eqn. (J^)], all one has to 
do is determine the Polyakov loop on each con- 
figuration separately (where it is a complex num- 
ber) , and then take the arithmetic average (which 
is supposed to be approximately real). The same 
holds true for the sectoral ensemble averages [i.e. 
the practical versions of the {L)\^\ showing up on 
the r.h.s. of (H)]. Hence, in the case of the small 
X simulation, instead of computing the sample 
average from the 3197 indexed configurations di- 
rectly (which gives L~ 0.495), one can also deter- 
mine (L)o and (L)i from the 3052, 145 configura- 
tions with [i^l =0, 1 (which gives 0.511 and 0.159) 
and then combine these sectoral averages with the 
weights po = 3052/3197, pi = 145/3197. In the 
run representing the intermediate LS-regime, the 
sectoral averages are (L)o = 0.361, (L}i — 0.228, 
(L)2 = 0.142, which together with the weights 
Po = 2343/3163, pi = 780/3163, pa = 40/3163 re- 
produce the ensemble average (L) =0.326. In the 
large x simulation, the corresponding numbers 
are (L)o = 0.0947, (L)i = 0.0861, (L)2 = 0.0777, 
(L)3 = 0.0003 and po = 924/2504, pi = 1253/2504, 
P2 = 293/2504, p3 = 34/2504, from which one also 
gets the ensemble average (L) = 0.0871. 

What one gains from this exercise is the in- 
sight that in the small x regime the physical ex- 
pectation value (L) is composed of sectoral av- 
erages {L)q, (L)i which prove very much incon- 
sistent, i.e. in the symmetry restoration regime 
a typical observable depends quite drastically on 
the topological charge of the background. In 
the intermediate LS-regime differences between 
neighboring topological sectors are much smaller, 
but the system is still far from showing over- 
all consistency among all {L)\^\. In the large 
X regime differences between neighboring topolo- 
gical sectors happen to be marginal, but ~ due 
to the highest topological sector with \v\ =3 - 
data indicate that even in the large LS-regime 
some limitations to the claimed |4J insensitivity 
on topology may persist. 



6. SUMMARY AND CONCLUSION 

The net outcome of the investigation presented 
here is that the predictions by Leutwyler and 
Smilga H] regarding the overall distribution of 
the topological charge v in the regimes a: <C 1 and 
a;3>l are well reproduced. Furthermore, the data 
ensure that the large x regime (where SSB of the 
axial flavour symmetry is pronounced) is special 
in several respects: On a formal level the situa- 
tion for a; :^ 1 is unique, as the functional determi- 
nant results only in a (mild) overall suppression 
of higher topological sectors w.r.t. lower ones, and 
the effective renormalization of /3 is uniform, i.e. 
independent of \v\. On a more phenomenological 
level, the analysis of the sectoral dependence of 
a typical observable has shown that the regime 
x ^ 1 is unique, since here physics proves (up to 
a limitation discussed above) independent of the 
topological charge of the background. Remarkab- 
ly, our results illustrate and extend those of 0|, 
even though in our x'^l and x^\ simulations 
the condition Mt^L ^ 1 in the LS-analysis has 
been reversed into A/jri^ 1, i.e. the pion would 
fit into the box. 
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